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(N- Abstract 
> ■ 

Q , We study the integrability of mappings obtained as reductions of the discrete Korteweg-de Vries 

■ (KdV) equation and of two copies of the discrete potential Korteweg-de Vries equation (pKdV) . We 

I show that the mappings corresponding to the discrete KdV equation, which can be derived from the 
latter, are completely integrable in the Liouville-Arnold sense. The mappings associated with two 

" copies of the pKdV equation are also shown to be integrable. 



^ ! 1 Introduction 

a- 

The problem of integrating differential equations goes back to the origins of calculus and its application 
to problems in classical mechanics. In the nineteenth century, the notion of complete integrability was 
provided with a solid theoretical foundation by Liouville, whose theorem gave sufficient conditions for a 

■ Hamiltonian system to be integrated by quadratures; yet only a few examples of integrable mechanical 
OO . systems (mostly with a small number of degrees of freedom) were known at the time. Poincare's 

I subsequent results on the non-integrability of the three-body problem seemed to indicate that the 

\^ • quest for other integrable systems was likely to be fruitless. Nevertheless, examples of integrable 

systems (and action-angle variables in particular) played an important role in the early development of 
quantum theory. 

■ The theory of integrable systems only began to expand rapidly in the latter part of the twentieth 
century, with the discovery of the remarkable properties of the Korteweg-de Vries (KdV) equation, 
together with a host of other nonlinear partial differential equations that were found to be amenable to 

X : the inverse scattering technique. As well as having exact pulse-like solutions (solitons) which undergo 

elastic collisions, such equations could be interpreted as infinite-dimensional Hamiltonian systems, with 
an infinite number of conserved quantities. Moreover, it was shown that these equations admitted 
particular reductions (e.g. to stationary solutions, or to travelling waves) which could be interpreted as 
integrable mechanical systems with finitely many degrees of freedom. The papers in the collection |37j 
provide a concise and self-contained survey of the theory of integrable ordinary and partial differential 
equations; for a more recent set of review articles, see |21j . 

In the last two decades or so, there has been a more gradual development of discrete integrable 
systems, in the form of both finite-dimensional maps [24:\ [5], discrete Painleve equations [26j, and 
partial difference equations defined on lattices or quad-graphs [22], [l] . Discrete integrable systems can be 
obtained directly by seeking discrete analogues of particular continuous soliton equations or Hamiltonian 
flows [30], but they also appear independently in solvable statistical mechanical or quantum models 
(see [18], for example). An important theoretical result for ordinary difference equations or maps is 
the fact that the Liouville-Arnold definition of integrability for Hamiltonian systems can be extended 
naturally to symplectic maps [191 HI El], so that an appropriate modification of Liouville's theorem 
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holds. For lattice equations with two or three independent variables there is much less theory available 
(especially from the Hamiltonian point of view) , and the full details of the known integrable examples 
are still being explored, but one way to gain understanding is through the analysis of particular families 
of reductions. 

By imposing a periodicity condition, integrable lattice equations can be reduced to ordinary differ- 
ence equations (or mappings) [Sj [251 13 Il6j- It is believed that the reduced maps obtained from an 
integrable lattice equation are completely integrable in the Liouville- Arnold sense. To prove that a map 
is integrable one needs to find a Poisson structure together with a sufficient number of functionally 
independent first integrals, and then show that these integrals commute with respect to the Poisson 
bracket. One complication that immediately arises is that the reduced maps naturally come in families 
of increasing dimension, and the number of first integrals grows with the dimension. Recently, some 
progress has been made with particular families of maps. For maps obtained as reductions of the equa- 
tions in the Adler-Bobenko-Suris (ABS) classification [Ij, and for reductions of the sine-Gordon and 
modified Korteweg-de Vries (mKdV) equations, first integrals were given in closed form by using the 
staircase method and the noncommutative Vieta expansion |3H I17j. In particular, the complete inte- 
grability of mappings obtained as reductions of the discrete sine-Gordon, mKdV and pKdV equations 
was studied in detail in [32l [33] . 

Given a map, the question arises as to whether it has a Poisson structure, and if so, how can one 
find it? In general, we do not know the answer. However, for some classes of maps one can assume that 
in coordinates xj the Poisson structure is in canonical or log-canonical form, i.e. the Poisson brackets 
have the form {xi,Xj} = Qij or {xi,Xj} = iljj XiXj respectively, where is a constant skew-symmetric 
matrix. This approach is valid for mappings obtained in the context of cluster algebras [121 M [IH] ; 
and also applies to reductions of the lattice potential KdV (pKdV), sine-Gordon and mKdV equations 
|33) . Another approach requires the existence of a Lagrangian for the reduced map: by using a discrete 
analogue of the Ostrogradsky transformation, as introduced in [4], one can rewrite the map in canonical 
coordinates; from there one can derive a Poisson structure in the original variables. 

Here we start by considering a well known integrable lattice equation, namely the discrete potential 
Korteweg-de Vries equation, also referred to as Hi in the ABS list [1], which is given by 

{Ue,m — Ue+i^rn+l){ue+l^m — U£,m+l) = Ij (1) 

where {£, m) G . Early results on this equation appear in [36] , [5] and [22] . As shown in [36] and [22] , 
the equation ([1]) leads to the continuous potential KdV equation, that is 

Ou_d^u^/dny^ (2) 



dt dx^ \dx J 

by performing a suitable continuum limit. Using the so-called three-leg form, as in [1], a Lagrangian 
can be obtained for the discrete pKdV equation ([T]). However, the associated Euler-Lagrange equation 
turns out to consist of two copies of the lattice pKdV equation, that is 

1 1 

Ug+l^m — Ue,m+1 + Ue-l^m — ue,m-i 1 = 0. (3) 

Ul,m — Ul+l,m+l Ul-l,m-l — Ui,m 

The latter equation is somewhat more general than ([!]) : every solution of ([T]) is a solution of ([3]) , but the 
converse statement does not hold. In this paper we shall be concerned with the more general equation 
§1, rather than 1^. 

Next, we introduce a new variable on the lattice, 

Ve,m '■= Ue,m — U£+l,m+l, 
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and immediately find that, whenever ug^m is a solution of v^^m satisfies 



1 



1 



(4) 



The latter equation is known as the lattice KdV equation. It is an integrable lattice equation, in the 
sense that it can be derived as the compatibility condition for an associated linear system, known as a 
Lax pair; this is discussed in section 4. 

In this paper we perform the {d — 1, — l)-reduction on the Lagrangian of the pKdV equation and 
derive the corresponding reduction of the Euler-Lagrange equation This means that we consider 
functions u = Ui^m on the lattice which have the following periodicity property under shifts: 



Such periodicity implies that u depends on the lattice variables i and m through the combination 



only; thus, with a slight abuse of notation, we write u = Un- Such a reduction can be understood as 
the discrete analogue of the travelling wave reduction of a partial differential equation: for a function 
u{x,t) satisfying a suitable partial differential equation such as ([2]), one considers solutions that are 
invariant under x ^ x + c6, t ^ t + 6 for all 6; such solutions depend on x, t through the combination 
z = x—ct only, corresponding to waves travelling with constant speed c. By analogy with the continuous 
case, where one obtains ordinary differential equations (with independent variable z) for the travelling 
waves, it is apparent that imposing the condition ([5|) in the discrete setting leads to ordinary difference 
equations (with independent variable n). Here we are concerned with the complete integrability of 
the ordinary difference equation obtained as the (d — 1, — l)-reduction of ([3]), which is equivalent to a 
birational map in dimension 2d, and the associated reduction of the lattice KdV equation which 
gives a map in dimension d. 

The paper is organized as follows. In section 2, we start with a discrete Lagrangian for the lattice 
equation ^ , corresponding to two copies of the lattice pKdV equation ([1]) . We then derive a symplectic 
structure for the reduced mappings obtained from the discrete Euler-Lagrange equation, by using a 
discrete analogue of the Ostrogradsky transformation. This provides a nondegenerate Poisson bracket 
for each of these maps. In section 3 we present a Poisson structure for the mapping obtained as the 
(d— 1, — l)-reduction of the lattice KdV equation for each d, which is induced from the bracket found 
in section 2. The Hirota bilinear form of each of these mappings is also given, in terms of tau- functions, 
whence (via a link with cluster algebras) we derive a second Poisson structure that is compatible with 
the first. In section 4 we present closed- form expressions for integrals of the reduced KdV map, fitting 
into the framework of the papers [31] and [33], which furnishes a direct proof of Liouville integrability. 
Another proof, based on the pair of compatible Poisson brackets, is also sketched. The construction of 
first integrals is done using the staircase method, the noncommutative Vieta expansion and a suitable 
splitting of Lax matrices. In section 5, for each d we present integrals for the mapping obtained as 
the {d — 1, — l)-reduction of the discrete Euler-Lagrange equation ([3]). Some integrals are derived from 
the integrals of the KdV map in the previous section, while additional commuting integrals for the 
reduction of ([3]) are found using a function periodic with period d; this is a d-integral, in the sense of 
[13] . The paper ends with some brief conclusions, followed by a proof that is relegated to an Appendix. 



U-l^m — Ui+d-l,m-l- 



(5) 



n = i -\- m{d — 1) 



3 



2 Nondegenerate Poisson structure from two copies of pKdV 

Henceforth, it is convenient to use ~ and ^ to denote shifts in i and m directions respectively, so that 
u = u = Ui^rn+ij 6tc. It is known that all equations of the form 

Q{u, u, u, n, a, /?) = (6) 

in the ABS list [Ij are equivalent (up to some transformations) to the existence of an equation of 
so-called three-leg type, that is 

P(u, u, u, u; a, f3) = <j){u, u, a) — (j){u, u, j3) — ip{u, u, a, (3) = 0, (7) 

for suitable functions (p, ip, where a and /3 are parameters; the latter leads to the derivation of a La- 
grangian for each of the equations ([6]). In particular, for the pKdV equation ([T]), which is a (parameter- 
free) equation of the form ([6|), we have (f){u,u) = u + u and 'ip{u,u) = l/{u — u), and using the three- leg 
form ([7]) leads to the Lagrangian 

£ = i(n + 5)2- l(u + S)2- login -n|. (8) 

The corresponding discrete Euler-Lagrange equation is the equation ([3]), which can be rewritten as 

J = J, with 3 = u — u —. (9) 

(u — u) 

The latter equation is more general than ([1]), which arises in the special case that J is identically zero. 

Now we set n = £ + m{d — 1), with u satisfying ([5]), and apply the {d — 1, — l)-reduction to the 
Lagrangian £ in ([8]), to find 

C {Un,Un+l, ■ ■ ■,Un+d) = ^{Un + Un+lf - ^{Un + Un+d-lf " log \Un - Un+d\- (10) 

The discrete action functional is S := Ylnei. ^ {un,Un+i, • • • , Un+d)- It yields the discrete Euler-Lagrange 
equation 

_ s;-^ dC {Un-r,Un+l-r, ■ ■ ■ ,Un+d-r) _ ST^ g^r £ _g Qj-v 



where £ denotes the shift operator and Cr = — g^^^^'"'""'*''' . Thus we obtain the ordinary difference 



dC{u„,U„+l,...,Un+d) 

equation 

Un+l - Un+d-1 + Un~l - Un~d+l \ = 0, (12) 

'^n '^n+d '^n—d 

which is precisely the {d — 1, — l)-reduction of The solutions of this equation are equivalent to the 
iterates of the 2(i-dimensional map 

(Un-d, Un~d+1, ■ ■ ■ , Un+d-l) {Un-d+l,Un-d+2, ■ ■ ■ , Un+d), (13) 

where Un+d is found from equation p^ . 

Given a Lagrangian of first order for a classical mechanical system, the Legendre transformation 
produces canonical symplectic coordinates on the phase space; the Ostrogradsky transformation is the 
analogue of this for Lagrangians of higher order [2j . In order to derive a nondegenerate Poisson bracket 
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for the 2(i-dimensional map, we use a discrete analogue of the Ostrogradsky transformation, as given 
in [1], which is a change of variables to canonical coordinates, 

{Un~d, ...,«„,... ,Un+d-l) {qi, ■ ■ ■ ,qd,Pl, ■ ■ ■ ,Pd), 

where 

d—i 

Qi = Un+i-1, Pi = S"^ ^ Cr+i- 

r=0 

Thus, from (fTOj) . we obtain 

qi = Un+i-i, i = l,...,d, 

1 1 

Pi = -Un+l + Un+d-1 H = Un-1 - Un+l-d + 



Pi = -Un+i-d - Un+i~l H ^ , i = 2, . . . , d - 1, 

^n-d+i-l ~ Un+i-l 

1 

Pd = • 

In terms of the canonical coordinates qj , pj the map (|12|) is rewritten as 

Qi H- Qi+i, i = l,...,d-l, Qd ^ qi - (Pi + q2 - Qd)'^ , 

Pi P2 + q2 + qi, Pi ^ Pi+1, i = 2,...,d-2, 

Pd-i ^ Pd-qi- qd, Pd ^ Pi + q2- qd- 

By the general results in [A], this map is symplectic with respect to the canonical symplectic form 
Yl'j=i '^Pj ^ dgj. Equivalently, it preserves the canonical Poisson brackets 



{pi,Pj} = {qi,qj} = 0, {pi,qj} = S: 



V 

In order to find the Poisson brackets for the coordinates we write them in terms of iqj,Pj)- For 
all < i < d — 1, we have Un+i = qi+i and 

Un-l = qd + ^IPd = [qd]Pd]-, 

Un-2 = qd-1 H ] — - — \ = [qd-i]Pd-i + qd-i + qd,Pd], 

Pd-l + qd+l + Un-l 
Un-i-1 = [qd-i\Pd-i + qd-i + Un-i] 

= [qd-i;pd-i + qd-i + qd-i+i,Pd-i+i + qd-i+i + qd-i+2, ■ ■ ■ ,Pd-i + qd-i + qd,Pd\, i < i < d - 2, 

Un-d = [qi;Pl + Un+l-d - Un-l] 

= [qi;pi -qd + q2- —,P2 + 1?2 + ga, • • • ,Pd-i + Qd-i + qd,Pd\, 
Pd 

where [ ; ] denotes a continued fraction. This yields the following result, proved in the Appendix. 



Theorem 1. The 2d- dimensional map given by I113\) with is a Poisson map with respect to the 
nondegenerate bracket 

{Un-d, Un-j-l} = 0, 0<j<d-l; (14) 
{Un-d,Un+j} = {-iy^^{Un-d " Unf ■ ■ ■ {Un-d+j " Un+jf , < j < d - 1] (15) 
{Un-d, Un+d-l} = (-l)'^('Un--d " Un)"^ ■ ■ ■ {Un-1 - Un+d-lf " {Un-d " Unf {Un-1 - Un+d-lf ■ (16) 
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3 Poisson structures and tau-functions for KdV maps 



As mentioned above, the discrete KdV equation can be derived from two copies of the pKdV equation. 
This suggests that we can find a Poisson structure for the {d — 1, — l)-reduction of the discrete KdV 
equation using the symplectic structure given in the previous section. It turns out that, in addition 
to the bracket induced from two copies of pKdV, each of the KdV maps has a second, independent 
Poisson bracket, which is obtained from a Hirota bihnear form in terms of tau-functions. The second 
bracket is constructed by making use of a connection with Somos recurrences and cluster algebras. 



3.1 First Poisson structure from pKdV 

Upon introducing the quantities 



we see that, whenever Un is a solution of equation p2|) . Vn satisfies a difference equation of order d, 
namely 

1 1 

Vn+d-l - Vn+1 \ = 0. (17) 

Vn+d "^n 

Alternatively, by starting from a solution of ([3]) with the periodicity property ([5|), we see that this yields 
a solution v = vi^m of with the same periodicity, and so (writing this as Vn, with the same abuse of 
notation as before) it is clear that equation pT]) is just the {d — 1, — l)-reduction of the discrete KdV 
equation. Equivalently, the ordinary difference equation (jl7p corresponds to the d-dimensional map 

(f : {vo,vi,...,Vd^i)^ ivi,V2,...,Vd-i,— — ). (18) 

The case d = 2 is trivial, so henceforth we consider d > 3. 

In the above, the suffix n has been dropped, taking a fixed d-tuple {vq, vi, . . . , Vd-i) in d-dimensional 
space. However, because the map (llSh is obtained from a recurrence of order d, all of the formulae are 
invariant under simultaneous shifts of all indices by an arbitrary amount n, i.e. Vj Vn+j for each 
j. For a fixed n, say n = 0, the formulae in Theorem [T] define a Poisson bracket in dimension 2d, 
which can be used to calculate the brackets between the quantities Vj = Uj_d ~ for j = 0, . . . ,d — 1. 
Remarkably, these brackets can be rewitten in terms of vj alone; in other words, these quantities form 
a Poisson subalgebra of dimension d. Hence this provides the first of two ways to endow (jlSp with a 
Poisson structure. 

Theorem 2. The d- dimensional map U8\) is a Poisson map with respect to the bracket { , }i defined 
by 

This bracket is degenerate, with one Casimir when d is odd, and two independent Casimirs when d is 
even. 

The above result follows immediately from Theorem[Tl apart from the statement about the Casimirs, 
which will be explained shortly. Here we first give a couple of examples for illustration. 

Example 3. When d = 3 the map ip given by (jlSp preserves the Poisson bracket 



{vq, Vi}i = -vlvl, {vi,V2}l = -vlvl, {vq,V2}i = {l+vl)vlv. 



2- 
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This bracket has rank two, with a Casimir C that is also a first integral for the map: 

C = vi , and ip*C = C. ^ 

Vq Vi V2 

Example 4. When d = A the map (llSp preserves the Poisson bracket specified by 

{^o,^i}i = -^^0^1, {vo,V2}i = vlvlvl, {vQ,Vi}i = {I - vlvl)vlvl, 

where all other brackets for < i,j < d — 1 are determined from skew-symmetry and the 

Poisson property of (/?. This is a bracket of rank two, having two independent Casimirs given by 

Ci=Vi-— C2 = V2-— ^, with (/3*Ci=C2, (p*C2=Ci. 

Vq V2 Vi Vs 

Remark 5. The Casimirs Cj in the preceding example are 2-integrals 113], meaning that they are 
preserved by two iterations of the map, i.e. {ip*)'^Cj = Cj for j = 1,2. The .symmetric functions 



JC = Ci+C2, JC' = CiC2 (20) 

provide two independent first integrals. 

In order to make the properties of the Poisson bracket { , }i more transparent, we introduce some 
new coordinates, the motivation for which should become clear from the Lax pairs in section 4. 

Lemma 6. For all d > 4, with respect to the coordinates 

go = -l/vo, gj = vj^i - 1/vj, j = 1, . . . , d - 1, (21) 

the first Poisson bracket for the KdV map U8\) is specified by the following relations for v = \ and 
<i < j <d-l: 

-1 ifj-i = l^ 
1 ifj-i = d-l, 
^/gl if i = '^,3 =d-l, 
otherwise. 



(22) 



When d is odd this bracket has the Casimir 



C = (?o + 9i + ---+5d-i + — , (23) 

90 

while for even d there are the two Casimirs 

Ci = 90 + 92 + ■ ■ ■ + gd-2, C2 = gi + g3 + ■ ■ ■ + gd~i + —■ (24) 

Remark 7. When v = Q the Poisson bracket \2S^) is just the first Poisson bracket for the dressing 
chain, as given by equation (13) in f35^ . 

In terms of the coordinates gj, the map (jlSp is rewritten as 

f- (go,gi,- • ■,9d-i) ^ (gi + —,92, ■ ■ -^gd-i, . ) ■ (25) 

V 50 1 + 9091 ) 

Note that for the special case d = 3, as in Example [3l part of the formula for the bracket ()22p requires 
a slight modification, namely 

{g\,g2\\ = -1 + i/fi-o- 
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3.2 Second Poisson structure from cluster algebras for tau- functions 

The discrete KdV equation was derived by Hirota in terms of tau-functions, via the Backlund 
transformation for the differential-difference KdV equation [H]. In Hirota's approach, the solution of 
the discrete KdV equation is given in terms of a tau- function as v = tt/{tt), and at the level of the 
{d — 1, — l)-reduction this becomes 

Vn = . (26) 

By direct substitution, it then follows that Vn is a solution of (jl7p provided that r„ satisfies the trilinear 
(degree three) recurrence relation 

Tn+2dTn+d-l Tn+1 = Tn+2d-l ^n+d+l Tn — 7"n+2d-l Tn+d-l Tn+2 + Tri,+2d-2 ^n+d+l Tn+l- (27) 

However, this relation can be further simplified upon dividing by Tn+iTn+di which gives the relation 

TnTn+d+l — Tn+2Tn+d-l Tn+d-lTn+2d — Tn+d+lTn+2d-2 

Oin = = = Oin+d-l- 

Tn+l Tn+d Tn+d 'Tn+2d-l 

This immediately yields relations that are bilinear (degree two) in r„. 



Proposition 8. The solutions of the equation are given in terms of a tau-function by I126\) . where 
Tn satisfies the bilinear recurrence relation 

'^n+d+l Tn = UnTn+dTn+l + T'n+d-l 'Tn+2 (28) 

of order d+\, with the coefficient a„ having period d — 1. 

Apart from the presence of the periodic coefficient a„, the bilinear relation (j28p has the form of a 
Somos-(d + 1) recurrence [29]. Such recurrence relations (with constant coefficients) are also referred 
to as three-term Gale-Robinson recurrences (after [llj and [28j respectively). 

Example 9. For d = 3 the equation ()28|) is a Somos-4 recurrence with coefficients of period 2: 

Due to the Laurent phenomenon the iterates of this recurrence are Laurent polynomials, i.e. poly- 
nomials in the initial values data and their reciprocals with integer coefficients; to be precise, 

Tn G 'L[ao,ai,T^^,Tf'^,T^^,T^^] 

for all n € Z. This means that integer sequences can be generated from a suitable choice of initial data 
and coefficients. For instance, with the initial values tq = ri = r2 = T3 = 1 and parameters ao = 1; 
ai = 2, the Somos-4 recurrence yields an integer sequence beginning with 

1, 1, 1, 1, 2, 5, 9, 61, 193, 1439, 13892, 121853, 1908467, 47166783, .... 

From the work of Fordy and Marsh [8], it is known that, at least in the case where the coefficients are 
constant, recurrences of Somos type can be generated from sequences of mutations in a cluster algebra. 
For the purposes of this paper, the main advantage of considering the cluster algebra is that it provides 
a natural presymplectic structure for the tau-functions. A presymplectic form that is compatible with 
cluster mutations was presented in pjj, and in [lOj it was explained how this presymplectic structure 
is preserved by the recurrences considered in [8]. 
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Cluster algebras are a new class of commutative algebras which were introduced in [B] . Rather than 
having a set of generators and relations that are given from the start, the generators of a cluster algebra 
are defined recursively by an iterative process known as cluster mutation. For a coefficient-free cluster 
algebra, one starts from an initial set of generators (the initial cluster) of fixed size, which here we take 
to be d + 1. If the initial cluster is denoted by (ti, . . . ,Tii+i), then for each index k one defines the 
mutation in the k direction to be the transformation that exchanges one of the variables to produce a 
new cluster (r(, . . . , t^_|_]^) given by 

d+l d+l 

r', = r„ j^k, r^, = nr^'" + n^'"'^"'"' (29) 

j=i i=i 

where the exponents in the exchange relation for come from an integer matrix B = (bij) known as 
the exchange matrix, and we have used the notation = max(6, 0). 

As well as cluster mutation, there is an associated operation of matrix mutation, which acts on the 
matrix B; the details of this are omitted here. Fordy and Marsh gave conditions under which skew- 
symmetric exchange matrices B have a cyclic symmetry (or periodicity) under mutation, and classified 
all such B with period 1 [8]. They also showed how this led to recurrence relations for cluster variables, 
by taking cyclic sequences of mutations. The requirement of periodicity puts conditions on the elements 
of the skew-symmetric matrix B, which (for a suitable labelling of indices) can be written as 

bi^d+i = h^i+i, i = l,...,d, (30) 

and 

bi+i,j+i = bij + - l<i,j <d. (31) 

The corresponding recurrence is defined by iteration of the map 

{ti,... ,Td,Td+l) ^ {T2,...,Td+l,T[) 

associated with the exchange relation (I29p for index k = 1, where the exponents are given by the entries 
bij in the first row of B. Moreover, given a recurrence relation of this type, the conditions (jSOp and 
()3ip allow the rest of the matrix B to be constructed from the exponents corresponding to the first 
row, and these conditions are also necessary and sufficient for a log-canonical presymplectic form u, as 
in (j32p below, to be preserved (see Lemma 2.3 in [lOj). In general this two-form is closed, but it may 
be degenerate. 

For the case at hand, the exponents appearing in the two monomials on the right hand side of (j28p 
specify the first row of the matrix B as (0, 1, —1, 0, . . . , 0, —1, 1), and the rest of this matrix is found 
by applying ([30l) and (fSTI) . Although the foregoing discussion was put in the context of coefficient- 
free cluster algebras, the presence of coefficients in front of these two monomials does not affect the 
behaviour of the corresponding log-canonical two-form under iteration. Thus we obtain 

Lemma 10. For all d > 5, the Somos-{d + 1) recurrence ^B^ preserves the presymplectic form 

oj = ^^dTi^ dTj , (32) 

given in terms of the entries of the associated skew- symmetric exchange matrix B = ibij) of size d+l, 
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where (with an asterisk denoting the omitted entries below the diagonal) 



( 1 





-1 



1 \ 

-1 





B 



2 




v 



Remark 11. In each of the special cases d = 3,4, the exchange matrix does not fit into the general 
pattern above. For the details of the case d = 3, when (28\) is a Somos-4 relation, see Example 2.11 in 
]10^ : and for d = 4 (Somos-5) see Example 1.1 in f^. 



The exchange matrices being considered here are all degenerate: for d odd, B has a two-dimensional 
kernel, while for d even, the kernel is three-dimensional. In the odd case, the kernel is associated with 
the action of a two-parameter group of scaling symmetries, namely 



while in the even case there is the additional symmetry 



(-1)" 



(33) 



(34) 



in the theory of tau- functions such scalings are known as gauge transformations. By Lemma 2.7 in |10j 
(see also section 6 therein), symplectic coordinates are obtained by taking a complete set of invariants 
for these scaling symmetries. With a suitable choice of coordinates, denoted below by yn, the associated 
symplectic maps can be written in the form of recurrences, like so: 



• Odd d: The quantities y„ 
equation 



Tn+2Tn/Tn+i invariant under ([33j) . and satisfy the difference 



Vn+d-l (yn+d-2 ' ' ' Vn+l) Vn = Oln yn+d-2 ' ' ' Vn+l + 1- 



(35) 



Even d: The quantities y„ = t„-(_3T„/(7:„+2T„+i) are invariant under both ([33]) and p^ . and 
satisfy the difference equation 



yn+d-2 ■ ■ ■ Vn = a-n yn+d-iUn+d-b 



Vn+l + 1- 



(36) 



Example 12. Both of the cases d = 3 and d = A lead to iteration of maps in the plane with coefficients 
that vary periodically, namely 

d='i: yn+2 Un+i Vn = Q-n Vn+l + 1, with an+2 = an, 

and 

d = A: yn+2 Un+i Vn = «n Un+i + 1, with an+3 = an- 

Remark 13. Due to the presence of the periodic coefficients, the latter maps are not of standard QRT 
type 1 24 J , but they reduce to symmetric QRT maps when a„ = constant. Maps in the plane of this more 
general type have recently been studied systematically by Roberts \2^ . 
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In general, the solutions of (|35p or (j36|) correspond to the iterates of a symplectic map, in dimension 
d — 1 or d — 2, respectively. To be more precise, rather than just iterating a single map, each iteration 
depends on the coefficient a„ which varies with a fixed period (as in Proposition [8|), but the same 
symplectic structure is preserved at each step. The appropriate symplectic form in the coordinates Un 
can be obtained from ([32|) . and by a direct calculation the associated nondegenerate bracket is found; 
this is presented as follows. 

Lemma 14. For d odd, each iteration of 1135]) preserves the nondegenerate Poisson bracket specified by 

{Vi, VA = i-iy-'^' ViVj, 0<i<j<d- 2. (37) 
For d even, the map defined by i36\) preserves the nondegenerate bracket in dimension d — 2 given by 

{yi,yi±i} = ±yiyi±i, (38) 

with all other brackets {yi, yj} for 0<i<j<d — 3 being zero. 

The bracket for the variables is the key to deriving a second Poisson structure for the KdV maps 
when d is odd. 



Theorem 15. In the case that d is odd, the map [TB^j preserves a second Poisson bracket, which is 
specified in terms of the coordinates [21\) by 



{9i,9ih = i-iy-'^^ 9i9j, 0<i<j<d- 1. (39) 



Proof. Substituting for from the formula (j26p and making use of the bilinear equation (j28p produces 
the identities 

^ 1 _ Tj{Tj+lTj+d-2 - Tj-lTj^d) _ aj-lTj 

SO that, by the definition of the symplectic coordinates yj, 

^ _aj_i^ J = 1, • • • - 1- 

yj-i 

A similar calculation in terms of tau-functions also yields 

90 = -yoyi ■ ■■yd-2- 

Noting that the coefficients play the role of constants with respect to the bracket (j37p . this immedi- 
ately implies that the induced brackets between the gi are 



{9i,9j} = oii-iaj-iy^\yj\{yi-i,yj-i} = (-1)^ ''^^gigj 

for l<i<j<(i — 1, which agrees with (j39p . By making use of the preceding formula for go in terms 
of yi, the brackets {5oi9j} follow in the same way. □ 

Remark 16. The bracket \3S^) is the same as the quadratic bracket for the dressing chain (for the case 
where all parameters j3i are zero in J3^). It has the Casimir 

C* = go gi ■ ■ ■ gd-i, (40) 
which is also given in terms of the coefficients of the bilinear equation \2^) by C* = — 17^=0 ■ 
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The case where d is even is shghtly more comphcated, because we do not have a direct way to derive 
a second Poisson bracket in terms of the coordinates gi (or equivalently Vi). The reason for this difficulty 
is that, from the tau-function expressions, although the quantities remain the same under the action 
of the two-parameter symmetry group (j33p . they are not invariant under the additional scaling (j34p . 
but rather they transform differently according to the parity of the index i: 

gi — > gi for even/odd i. (41) 

In order to obtain a fully invariant set of variables, we introduce a projection vr from dimension d to 
dimension d — 1: 

: fi=gigi+i, i = 0,...,d-2. 

By regarding the new coordinates /, as functions of the gj, we get an induced map (/?' in dimension 
d — 1, which is compatible with in the sense that (p ■ n = n ■ ip' ; this has the form 

V>' ■■ (/O, /l, . . . , fd-2) ^ (fl{l + fo')J2, fd-2, ,^ t-^'" 7''"% ) ■ (42) 

V {^ + Jo )J1J3 - ■ ■ Jd-3/ 

In terms of tau-functions, the quantities fi are invariant under the action of the full three-parameter 
group of gauge transformations, which means that they can be expressed as functions of the invariant 
symplectic coordinates yn, and hence we can derive a Poisson bracket for them. 

Theorem 17. In the case that d is even, the map preserves a Poisson bracket, which is specified 
in terms of the coordinates fi by 

{f^J^±l}2 = ±f^f^±l, (43) 

with all other brackets {fi, fj}2 for 0<i<j<d — 2 being zero. 
Proof. Following the proof of Theorem [T5l we have 

Ui-iaiTiTi+i 

9i9i+i — , 

Ti-l'Ti+2 

whence, in terms of the symplectic coordinates for (j36p . we have 

ai-iai . 
fi = , i = l,...,d-2. 

Vi-i 

A similar calculation yields a slightly different formula for index i = 0: 

fo = aoyiys ■ ■ ■ yds, 

and then the Poisson brackets (j38p between the yj directly imply that the brackets (j43p hold between 
the fi. □ 

Remark 18. The bracket ^3\ ) has the Casimir 

C* = fof2---fd-2, (44) 

which is also given in terms of the coefficients of the bilinear equation \28\) by C* = 11^=0 ^j- 

It may be unclear why there is the suffix 2 on the bracket in (j43p . since we have not yet provided 
another Poisson bracket for the map (I42p . However, as will be explained in the next section, when d 
is even the quantities fi form a Poisson subalgebra for the bracket {, }i of Theorem [2j As will also 
be explained, the brackets { , }i and { , }2 are compatible (in the sense that any linear combination of 
them is also a Poisson bracket), which means that a standard bi-Hamiltonian argument can be used to 
show Liouville integrability of the maps for either odd or even d. 
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4 First integrals and integrability of the KdV maps 



The purpose of this section is to prove the fohowing result. 

Theorem 19. For each d > 3, the map il8\) is completely integrable in the Liouville sense. 

Recah that, for a Poisson map (/? in dimension d, Liouville integrability means that there should be 
k Casimirs invariant under the map, so that the {d — /i:)-dimensional symplectic leaves of the Poisson 
bracket are preserved by ip, plus an additional ^{d — k) independent first integrals that are in involution 
with respect to the bracket. For the particular map in question, there is always the Poisson bracket 
{ , }i, and from Lemma E] this has either one or two Casimirs, with symplectic leaves of dimension d—\ 
or d — 2, for odd/even d respectively; hence an additional [^(d — 1)J first integrals are required in this 
case. 

For the lowest values li = 3, 4, it is straightforward to verify complete integrability, since in those 
cases only one extra first integral is required, apart from the Casimirs of {, }i. When d = 3, as in 
Example [3l the Casimir C is preserved by ip; the quantity C* in (j40p . which is the Casimir of the second 
bracket, is also a first integral. Similarly, for d = 4, the two Casimirs in Example H] provide the two first 
integrals ([20]) . which are themselves Casimirs; the quantity C* = /0/2 = 90919293, as in dM]), provides 
the extra first integral. 



Example 20. For d = 3, each of the first integrals C, C* of the map (jlSp define surfaces in three 
dimensions, given by 

XY'^Z - C XYZ - XY -YZ - ZX = 0, XY'^Z + C* XYZ - XY -YZ + 1 = 0, 

respectively, in terms of coordinates {X,Y,Z) = {vn,Vn+i,Vn+2)- These two surfaces intersect in a 
curve of genus one, found explicitly by eliminating the variable Z above; this yields the biquadratic 

(C + C*) X'^Y^ + {X + Y){XY - I) - C XY = 0. (45) 

The embedding of such a curve in three dimensions can be seen from the orbit plotted in Figure 1. <i) 

For higher values of d, it is not so obvious how to proceed, but the correct number of additional 
first integrals can be obtained by constructing a Lax pair for the map, as we now describe. 

4.1 Lax pairs and monodromy 

The discrete KdV equation dH) is known to be integrable in the sense that it arises as the compatibility 
condition for a pair of linear equations. In [23], a scalar Lax pair is given as follows: 

(f> = v(f> + X(f>, (j) = (j) -\ — (j). 

V 



Upon introducing the vector <^ := ( , the latter pair of scalar equations leads to a Lax pair in 
matrix form, namely 

$ = L$, 8 = M$, with L = ~ ^ , M = (46) 

Equation is equivalent to the compatibility condition for the linear system (|46p . that is 

LM = ML. (47) 
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Figure 1: The orbit of the map ([T^ for d = 3 with initial data {vo,vi,V2) = (1, 5, |). 



Note that the matrix L is associated with shifts in the i (horizontal) direction, and M with shifts in 
the m (vertical) direction. 

First integrals of each KdV map p8|) can be found by using the staircase method [5l [23 [151 [TB] . For 
the {d — 1, — l)-reduction, a staircase on the lattice is built from paths consisting of d — 1 horizontal 
steps and one vertical step. Taking an ordered product of Lax matrices along the staircase yields the 
monodromy matrix 

Cn = M„ ^ L„+(i_2 • • • L„, (48) 

corresponding to d — I steps to the right {i ^ i + 1) and one step down (m — > m — 1), where L — > L„ 
and M M„ under the reduction. As a consequence of (flTl) . the identity L„+rf_iM„ = M„+iL„ 
holds, which implies that satisfies the discrete Lax equation 

Cn+l- (49) 

The latter holds if and only if Vn satisfies the difference equation (I17p . Since (I49p means that the 
spectrum of Cn is invariant under the shift n — )• n+ 1, first integrals for the KdV map can be constructed 
from the trace of the monodromy matrix (or powers thereof), which can be expanded in the spectral 
parameter A. 

For convenience, we conjugate Cn in (US]) by M„ and multiply by an overall factor of A — 1, which 
(upon setting n = 0) gives a modified monodromy matrix 

^^ = LLr--LlLj, (50) 

where 

Ll = U^i=(^^l i = l,...,d-l, but lJ=(^]° i) with iy = l. 
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Thus we see the origin of the coordinates gi introduced in Lemma El they are the (1, 1) entries in the 
Lax matrices that make up the monodromy matrix. We have already mentioned a connection with 
the dressing chain at the level of the Poisson brackets, but it can be seen more directly here: up to 
taking an inverse and inserting a spectral parameter, when 1^ = the expression (j50p reduces to the 
monodromy matrix for the dressing chain found in [TO]. The parameter can be regarded as introducing 
inhomogeneity into the chain at i = (and for i/ ^ it is always possible to rescale so that = 1). 

With the introduction of another spectral parameter n, we consider the characteristic polynomial 
of C\ which defines the spectral curve 

X(A, fi) := det(/:t - ^ 1) = ^2 _ p(^) ^ + A)(-A)'^-i = 0, (51) 

with V{X) = tr C"^ . This curve in the (A,/i) plane is invariant under the KdV map ()18p . and the trace 
of £^ provides the polynomial V whose non-trivial coefficients are first integrals of the map. 

Example 21. When d = 3 the spectral curve is 

^2 - (C* + CX)fi + {u- A)A2 = 0, 

which has genus one, and (for z/ = 1) is isomorphic to the biquadratic curve (|45p corresponding to the 
intersection of the level sets of C and C*. <C> 

Example 22. When d = 4 the spectral curve also has genus one, being given by 

fi^ - (lo + /i A + 2A2) ^l-ily- A)A3 = 0, 
with coefficients expressed in terms of gi by 

Iq = C* = gogmgs, h = fC' -i^ = gogi+ gm +9293 +9390 + 1^ —■ 

90 

In the next subsection, we will show how, for each d, expanding the trace of the monodromy matrix 
in powers of A gives the polynomial V of degree \_d/2\. This implies that the hyperelliptic curve 
x(A,/i) = as in (f5T]l has genus [|(d — 1)J. Thus we expect that the (real, compact) Liouville tori for 
the KdV map should be identified with (a real component of) the Jacobian variety of this curve, since 
their dimensions coincide, and the map should correspond to a translation on the Jacobian. 



4.2 Direct proof from the Vieta expansion 

In order to calculate the trace of the monodromy matrix explicitly, we split the Lax matrices as 

L| = a(o + °)=AX + Y„ ^ = l,...d-l, (52) 

and similarly for Lq. This splitting fits into the framework of [31j, with the coefficient of A being the 
same nilpotent matrix X in each case, and allows the application of the so-called Vieta expansion. 
Recall that, in the noncommutative setting, the Vieta expansion is given by the formula 



J](AX, + Yi) := (AX6 + Y6)...(AX„ + YJ = ^ A^Z: 

i=a 

where 



b-a+l 

a,b 
r 1 

r=0 



Z"'*" — ^ Yf,Yf,_i • • • Yj,,+iXj^Yj,,_i • • • Yj^+iXj^Yj^_i • • • Ya. (53) 

a<ii<i2<-<ii<f) 
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In the case at hand, we have Xj = X for ah i, and by writing the monodromy matrix in the form (|53p . 
we can use Lemma 8 in |31j, so that the trace is expanded as 



K2J 

ViX) = trC^ = ^ /r■A^ (54) 

r=0 

where the first integrals of the KdV map are given by the closed-form expression 

T ,T,l,rf-3 , ,T,l.f^-4 , I I ^ \ ,T,2,rf-3 , ,T,2,(i-4 , ,t,14-3 /'c-r\ 

Ir = '^r-1 +90'^rU ^ i ^^^^^ ^ ^ ) + +909d-l'^r ) (55) 

with 

i=a \a<ii<ii+l<i2...<ir<b j=l ^*J^*J+^y 

The explicit form of these integrals mean that it is possible to give a direct verification that they are 
in involution, which is almost identical to the proof of Theorem 13 in |33] . 

Proposition 23. The first integrals ( 155|) of the KdV map Poisson commute with respect to the bracket 
{,}i in 

Proof. Writing the Poisson structure in terms of the coordinates gi as in Lemma [6l we see that (up to 
an overall factor of —1) the brackets {gi,gj}i in (I22p for 1 < i,j < d — 2 are identical to the brackets 
between the coordinates Cj given by equation (46) in [33j. The polynomial functions given by ()56p are 
the same as in the latter reference, and the particular functions ^"'^ that appear in the formula (j55p 
only depend on gi for 1 < i < d — 2. This implies that the brackets between these functions are the 
same as those given in Lemma 11 and Corollary 12 in ^33j, whence it is straightforward to verify that 
{Ir,Is}i = for < r, s < [d/2\ . □ 

In order to complete the proof of Theorem [191 it only remains to check that there are sufficiently 
many independent integrals. 

For odd d, the leading = 0) part of the first integral Ir is a cyclically symmetric, homogeneous 
function of degree d — 2r in the gi, being identical to the independent integrals of the dressing chain in 
[35j (for parameters /Jj = 0), with terms of each odd degree appearing for r = 0, . . . , Hence these 
functions are also independent in the case i' ^ 0, corresponding to the addition of a rational term with 
go in the denominator, which is linear in z/ and of degree d — 2r — 2 in the gi. This means that there 
are degV = (d — l)/2 + 1 independent integrals, of which the last one is I(d-i)/2 = C, the Casimir of 
the bracket { , }i, as given in ([23l) . 

When d is even, the leading (u = 0) part of Ir has the same structure, with terms of each even 
degree d — 2r appearing for r = 0, . . . , | — 1, but for all u the last one is trivial: 1^/2 = 2. Thus the 
coefficients of V provide only ^ independent first integrals. The last non-trivial coefficient is a Casimir 
of { , }i, given by Id/2-1 = C1C2 — i^-, in terms of the two Casimirs in (|24]l . and the Casimir Ci + C2 
provides one extra first integral, as required. 

In the next subsection, we briefiy outline another proof of Theorem [191 based on the bi-Hamiltonian 
structure obtained from the two Poisson brackets { , }i and { , }2- 

4.3 Proof via the bi-Hamiltonian structure 

Rather than making a direct calculation of brackets between the first integrals (j55p . the two Poisson 
brackets can be used to show complete integrability, by employing a standard method for bi-Hamiltonian 
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systems PO]. This approach apphes immediately to the coordinates gi when d is odd and the brackets 
are given by (j22p and (j39p : but minor modifications are required to apply it to the coordinates fi when 
d is even. 

For d odd, the first observation to make is that, just as in the case of the dressing chain (z^ = 0), 
the brackets { , }i and { , }2 are compatible with each other, in the sense that their sum, and hence 
any linear combination of them, is also a Poisson bracket. Thus these two compatible brackets form a 
bi-Hamiltonian structure for the KdV map 93. A pencil of Poisson brackets is defined by the bivector 
field P2 — A Pi, where Pj denotes the bivector corresponding to { , }j for j = 1, 2. (In fact, since it defines 
a bracket for all A and z^, this gives three compatible Poisson structures.) Then, by a minor adaptation 
of Theorem 4.8 in [lOj (which, up to rescaling by a factor of 2, is the special case A = —1 and = 0) 
we see that the trace of the monodromy matrix is a Casimir of the Poisson pencil, or in other words 

(P2-APi)jdP(A) =0. (57) 

Expanding this identity in powers of A yields a finite Lenard-Magri chain, starting with Iq = C*, the 
Casimir of the bracket { , }2 (as in Remark 1 16p . and ending with I[d-i)/2 = C, the Casimir of the bracket 
{ , }i (as in Lemma E]): 

P2jd/o = 0, P2 jd/i = Pud/o, ••• P2 Jd/r =PijdI^_i, ... Pijd/(d_i)/2 = 0. 

It then follows, by a standard inductive argument, that {Ir,Is}i = = {Ir, Is}2 for < r,s < (d—l) /2. 
Hence we see that Proposition [23] is a consequence of the bi-Hamiltonian structure in this case. 

In the case where d is even, a slightly more indirect argument is necessary, making the projection 
TT and working with the map ip' in dimension d — 1, as given by (j42p . The main ingredient required is 
the expression for the relations between the coordinates fi for i = 0, . . . ,d — 2 with respect to the first 
bracket. The case d = 4 is special, so we do this example first before summarizing the general case. 

Example 24. For d = 4, we use Lemma [6] to calculate 

{/o,/i}i = {9ogi,9i92}i = {go,9i}i9ig2 + gii{9o, 92)191 + 9o{9i,92}i) = -9192 - 9o9i = -/o/i, 
and similar calculations show that 

{/o, /2}i = fl - + 2^ T"' {/i' hh = - fl - f2 + v ^, 
h Jo Jo 

which implies that /o,/i,/2 generate a three-dimensional Poisson subalgebra for the bracket { , }i. It 
follows that (I42p is a Poisson map (in three dimensions) with respect to the restriction of this bracket 
to the subalgebra. Moreover, the restricted bracket { , }i for the fi is compatible with the bracket { , }2 
given by ()43p with d = A. The coefficients of the corresponding spectral curve, as in Example [22] can 
be written as functions of the fi, i.e. 

/0 = /0/2, Il=f0 + fl + f2+^-4^ + l^fy. 

Jl Jo 

and these generate the short Lenard-Magri chain 

{•,/0}2 = 0, {.,/i}2 = {.,/o}l, {•,/l}l=0. 
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Theorem 25. For all even d > 6, the first Poisson structure for the KdV map \18\) reduces to a bracket 
for the coordinates fi = gi Qi+i, as specified by the following relations for v = 1 and 0<i<j<d — 2: 



— fi — fi+1 

_ fifi+2 
fi+1 

/l/3---/d-3 
f2fi---fd-A 

y /l /3---/d-3 

/o/2---/d-4 





/o 



if j 



if 3 



1, 



d-2, 



(58) 



= l,i = d - 2, 



otherwise. 



This Poisson bracket is preserved by the map 

The rest of the argument proceeds as for d odd: the brackets { , }i and { , }2 (as given in ([58]) and 
(I43p respectively) are compatible with each other, so they provide a bi-Hamiltonian structure for ip' in 
dimension d — 1. Letting for j = 1,2 denote the corresponding Poisson bivector fields, the analogue 
of ()57p can then be verified: 

(P'2-AP'i)jd7'(A)=0. 

It should be noted that the latter identity is well-defined in the {d — l)-dimensional space with coor- 
dinates fi'. the trace of the monodromy matrix is invariant under the scaling symmetry (I4ip . hence all 
of the first integrals Ir can be written as functions of the variables fi. In this case, the Lenard-Magri 
chain begins with Iq = C* , the Casimir of { , }2 (as in Remark llSp . and ends with a Casimir of the 
bracket { , }i, namely Id/2-1 — C1C2 — i' (which is well-defined in terms of fi). Thus it follows that the 
integrals Ir are in involution with respect to both brackets for the fi, and this result extends to the 
bracket { , }i when it is lifted to d dimensions (in terms of gi, or equivalently Vi). 



5 Additional integrals for two copies of the pKdV equation 

In this section, we go back to the {d — 1, —1) reduction of the equation ([3]). This reduction yields the 
difference equation (jl2p . and in section 2 we showed how writing this as the discrete Euler-Lagrange 
equation (fTT]l led to a symplectic structure for the corresponding 2(i-dimensional map (fT3l) . It is easy to 
see that the first integrals (j55p of the KdV map also provide first integrals for equation (jl2p by writing 
Vi = Ui-d — «i for i = 0, . . . , d — 1. However, the total number of independent integrals for (fTSl) is only 
\_d/2\ + 1, which is not enough for complete integrability of the 2(i-dimensional symplectic map (jl3p . 
Here we will show how to construct sufficiently many additional integrals, leading to a proof of the 
following result. 

Theorem 26. For all d > 3, the 2d- dimensional map given by ( fjgj) with U^) is completely integrable 
in the Liouville sense. 

The proof of the above relies on the very simple observation that equation ()12p can be rewritten as 

1 1 

Un+l — Un+d-l — — ""n-d+l " Un-1 — , 

from which we see that 

Un-d+l - Un-1 ^3 = Jn, with Jn+d = Jn- (59) 
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Thus the function J„ is a d-integral (in the sense of [T3]) for the map for any n it can be viewed 

as a function on the 2d-dimensional phase space with coordinates {u-d, ■ ■ ■ ,Ud-i), and under shifting 
n it is periodic with period d. This imphes that any cyclically symmetric function of Jq, Ji, . . . , Jd-i is 
a first integral for equation (jl2p . As we shall see, for the corresponding map (jl3p this has the further 
consequence that it is superintegrable, in the sense that it has more than the number of independent 
integrals required for Liouville's theorem. 

Remark 27. For the (d— 1, —1) -reduction, the foregoing observation is a direct consequence of the fact 
that (0j can be rewritten as An analogous observation applies to the {si, —S2) -reduction of 
where n = + sim, and in ^ one has J — )• J„ with J„_|_j,j_|_s2 = Jn- 



5.1 Construction of integrals in involution 

In order to construct additional integrals, we need to calculate the Poisson brackets between the Jj, 
as well as their brackets with the Vj] together these generate a Poisson subalgebra for the bracket in 
Theorem [H as described by the next two lemmata. 

Lemma 28. With respect to the nondegenerate Poisson bracket specified by il5\) and U6\) . the 

quantities Ji defined in Ii59\) and Vj = Uj-d — Uj Poisson commute, i.e. 



{Ji,Vj} = 0, 0<i,j<d-l. (60) 

Proof. For < j < d — 1 we have 

{Jo,Vj} = {U-d+l - U_i - VQ^,Vj} = {U^d+l - U-l,Uj^d - Uj} + VQ'^{vo,Vj} 
= -{u.d+l,Uj} + Vq^ {-ly Ui=0 = 0, 

where we have used (fT9]l as well as Theorem [H Similar calculations show that { Jo, fo} = = {Jo,Vd-i}, 
and, since the bracket is preserved by the map (jl3p . the vanishing of all the other brackets {Ji,Vj} 
follows by shifting indices. □ 

Lemma 29. The Poisson brackets between the d-integrals Jn are specified by 

{1^ ij j — i = I 
-1, ifj-i = d-l, (61) 
0, otherwise, 

for < i < j < d - 1. 

Proof. Once again, from the behaviour under shifting indices, it is enough to verify the brackets for 
i = and I < j < d — 1. Expanding the left hand side of (j6ip . and using Lemma [281 '^6 obtain 

{Jo, Jj} = {Jo,Uj-d+l - Uj-i - Vj^} = {u^d+l - U-i - {U-d - Uo)'^ ,Uj^d+l - Uj-l] 

= -{n_d+i,Uj_i} + (n_d - uo)~^ {{u-d,Uj^d+i] - {u-d,Uj-i}^. 

Upon substituting with the non-zero brackets for the coordinates ui, as in Theorem [H the above 
expression vanishes for 2 < j < d — 2, and is equal to ±1 for j = 1, d — 1 respectively. □ 

The brackets ([6T]) mean that a suitable set of quadratic functions of the Jj give additional integrals 
for the equation (fT2]) . 
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Proposition 30. The functions 

Ts = Y,J^ Ji+s (62) 

1=0 

provide \d/2\ + 1 independent first integrals for the map { corresponding to the {d — 1, —\)-reduction 
of two copies of pKdV equation. These integrals are in involution with each other, and Poisson commute 
with the first integrals of the KdV map 17^) . 



Proof. With indices read mod d, the quantities Tg are cychcally symmetric functions of the J, ; in other 
words they are invariant under the cychc permutation (Jq, . . . , Jd-i) ^ (Ji, ■ ■ ■ , Jd-i, Jo), which means 
that they are first integrals for (jl3p . From the periodicity of the d-integrals Ji it is clear that Tg^a = Tg, 
and also Td-s = Tg. Taking s = 0, . . . , [d/2\ yields [d/2\ + 1 independent functions of Jq, . . . , Jd-i, and 
since the Ji are themselves independent functions of u_d, . . . , Ud-i, this implies functional independence 
of this number of the quantities ()62p . The fact that these quantities Poisson commute with each other 
is a consequence of Lemma [29l and the fact that 

oJi 

which implies that 

^ , ^ fmm_ _ _dT^dTA _ . 



as required. Also, by Lemma \28\ each Tg Poisson commutes with any function of the Vi, so with the 
first integrals of (jlSp in particular. □ 



Before describing the general case, we now demonstrate complete integrability for the simplest 
examples. 



Example 31. Starting with d = 3, equation (|12p yields the 6-dimensional symplectic map 

(n^3, n_2, ...,U2)'-^ (^n„2, ... ,U2,F (n_3, m_2, . . . , 1x2) ^ , (63) 



where 

F(li_3,U_2, ...,U2)=Uo 



Ui + U-i - U2- U-2 + l/(^t-3 - ^^o) ' 

Two integrals of this map are given in terms of vq = m_3 — wq, vi = U-2 — ui, V2 = ti-i — U2 by 
Io = C =-vi-\ 1 , Ii=C = vi 

Vo V2 V0V1V2 Vq Vi V2 

Apart from these, there is the pair of integrals 

To = Jq + Ji + J21 ^1 = •^•^1 + J1J2 + J2Jo-, 

which are written as symmetric functions of the 3-integrals Ji = Ui-2 — Ui-i — {ui-^^ — Ui)~^ , i = 0, 1, 2. 
The quantities Iq, Ii,Tq,Ti Poisson commute with each other, but they cannot all be independent. 
Indeed, the quantities vq,vi,V2, Jq, Ji, J2 are not themselves independent functions of Uj: they are 
connected by the relation 

h=vi = Jo + Ji + J2, 

Vq Vi V2 
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which imphes that the first integrals satisfy 

lf = To + 2Ti. 

Subject to the latter relation, the Liouville integrability of the map ()63p follows by taking any three 
independent first integrals in involution {Iq,Ii,To, for instance). The existence of an additional inde- 
pendent first integral, namely 

•^o = J0J1J2 

(another symmetric function of Jo, Ji, J2), means that the map is superintegrable; but this integral does 
not Poisson commute with Tq or Ti. <0 



Example 32. For d = 4, equation (|12p gives the 8-dimensional map 

(U_4, U_3, ...,U3)^ (u_3, ...,U3,G (u_4, U_3, . . . , ^3) ^ , (64) 

where ^ 

G {u_4, u_3, . . . , Us) = no ■ —77 -. 

Ui + U-l -U3 - U-3 + l/(n_4 - Uo) 

This map is symplectic with respect to the nondegenerate Poisson structure given in Theorem [TJ There 
are three independent integrals that come from the KdV map, given by 

lo = - 1 ) (vi - - ) (v2 --], h= C1C2 - 1, /C = Cl + C2, 



which are all written in terms of Vi = Ui-4 — Ui for i = 0, 1, 2, 3, using 

11^ 11 

Cl=Vl , C2 = V2 . 

Vq V2 Vi V3 

The first two of these integrals (/o and Ii ) are coefficients of the spectral curve in Example [22] (with 
1/ = 1), while the third is not. As well as these, there are three independent cyclically symmetric 
quadratic functions of the 4-integrals Ji = Ui-3 — Ui-i — (uj_4 — Ui)~^, i = 0, 1, 2, 3, namely 

To = Jq + Jf + J|, Ti = JqJi + J1J2 + J2J3 + T2 = 2(Jo J2 + J1J3), 



which are also first integrals of (|64p . There are two relations between the quantities vq, . . . ,V3 and 
Jo , . . . , J3 , as can seen by noting the identities 

Cl = Jq + J2, C2 = Ji + J3; 

thus the aforementioned first integrals are related by 

= ro + 2Ti + T2, /i + l = Ti. 



Hence complete integrability of the map (I64p follows from the existence of 4 independent integrals in 
involution, i.e. Iq, Ii,IC,Tq. The map is also superintegrable, due to the presence of a fifth independent 
first integral, given by another symmetric function of the Jj, that is 

5*0 = JqJi J2J3- ^ 
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As we now briefly explain, the general case follows the pattern of one of the preceding two examples 
very closely, according to whether d is odd or even. 

When d is odd, the spectral curve ()5ip has (d+l) /2 non-trivial coefficients, which are the quantities 
Ir appearing in ([5l|) . There are also {d + l)/2 independent functions Tg, as in (f62|) . but the identity 
^{d-i)/2 = Jq + Ji + ■ ■ ■ + Jd-1 implies that these two sets of functions are related by 

d-l 

^{d-l)/2 =Z]^^- 
s=0 

Hence there are precisely d independent functions in involution, as required for Theorem 1261 

In the case that d is even, the non-trivial coefficients of the spectral curve are supplemented by 
the additional integral IC = Ci +C2, providing a total of {d + 2)/2 independent functions, and there are 
the same number of independent functions of the form (I62p . but now the pair of identities 

Ci = ^ Ji, C2 = ^ Ji 

i even, 0<j<d— 2 jodd,l<i<d— 1 

together imply that the first integrals satisfy the two relations 

d-i 

= ^ Ts, I{d-2)/2 + 1 = ^ Ts, 

s=0 sodd,l<s<d/2 

SO once again there are d independent integrals, as required. 

One can also construct extra first integrals Sj for j = 0, . . . , [{d — l)/2j — 1, by taking additional 
independent cyclically symmetric functions of the Ji. This means that the map (|13p is superintegrable 
for all d. 



5.2 Difference equations with periodic coefficients 

To conclude this section, we look at equation (j59p in a different way, and show how it is related to other 
difference equations with periodic coefficients. We begin by revisiting the previous two examples. 

Example 33. For the map (f63l) . the introduction of the variables Xn = Un-3 — Un-2 yields the following 
difference equation of second order: 

Xn+2 + Xn = ^^—f J„+3 = Jn- (65) 

Xn+1 •Jn 

Iteration of the latter preserves the canonical symplectic form dxQ A dxi in the {xq,xi) plane. Apart 
from the coefficient J„, there is another periodic quantity associated with this equation, namely the 
3-integral Hn given by 

= X^Xn+l + x'^^iXn + {Jn+2Jn — Jn+lJn — ^)Xn + {Jn+2Jn — Jn+lJn+2 — ^)Xn+l 
JnX^ Jn+2Xj^^l + {Jn+l Jn Jn+2^XnXn-\-l 

(with Hn-\-3 = Hn)- This is related to some of the first integrals for d = 3 by the identity 

Hn — Jn+l = Iq — Sq, 

where Sq = J0J1J2 as in Example [3T1 ^ 



22 



Example 34. By introducing Wn = Un-i — Un-2 in we obtain the second order equation 

Wn+2 +Wn = ^- ^, Jn+A = Jn- (66) 

Wn+1 — Jn 

Each iteration of this difference equation preserves the canonical symplectic structure dzuo A dwi . Aside 
from the coefficient the equation (|66p has a 4-integral i.e. Hn+A = -f^n where 

+ {JnJn+2Jn+?, — JnJn+lJn+2 — Jn+2)Wn + {JnJn+lJn+3 — Jn+lJn+2Jn+3 — Jn+l) Wn+1 
— JnJn+2'w'^ — Jn+lJn+S'Wn+l + {Jn+3 — Jn+l) {Jn — Jn+2) — ^^WnWn+l- 

One can check that Hn+i — Hn = Jn+2{Jn+3 — Jn+i), and Hn is related to the first integrals in Example 
[32] by the formula 

Hn — Jn+l Jn+2 = Iq — Sq. <0' 

Remark 35. Both equations i65\} and i66\) are of the type considered recently by Roberts \21^ : their 
orbits move periodically through a sequence of biquadratic curves, defined by the quantities Hn, and they 
reduce to symmetric QRT maps when Jn = constant. 

In general, equation ()59p can be rewritten in terms of the variables Xn = Un-d — Un-d+i, to obtain 
a difference equation of order d — 1, that is 

Xn + Xn+1 + . . . + Xn+d-1 = : , ^ , (67) 

Xn+1 + Xn+2 + . . . + Xn+d-2 - Jn 

where J„ is periodic with period d. Note that both sides of the above equation are equal to Vn, the 
discrete KdV variable, as in ()17p . This equation can be seen as a higher order analogue of the McMillan 
map, with periodic coefficients. First integrals for equation ()67p can be obtained from the integrals Ir 
of the KdV map (jlSp by rewriting all the variables Vi (or gi) in terms of x„ and J„. To be precise, one 
can verify that gi = + a;, + Jj for 1 < i < d — 1, and vo = —gQ^ is given by the right hand side 
of (|67p for n = 0. In particular, the explicit formula for Iq is 

d-2 / /d-2 \ /d-3 \ \ 

Iq = Wixi-i +Xi + Ji){ y^Xi- Jq \ y^Xi- Jd-i - 1 • 

When d is even, one can reduce ([67|) to a difference equation of order d—2, by setting Wn = Xn+Xn+i- 
Formulae for first integrals in terms of the variables Wn follow directly from integrals given in terms of 
Xn and Jn- 



6 Conclusions 

We have demonstrated that the {d — 1, — l)-reduction of the discrete KdV equation ([3]) is a completely 
integrable map in the Liouville sense. There are two different Poisson structures for this map: one was 
obtained by starting from the related equation ^ and its associated Lagrangian; the other arose by 
using tau- functions and a connection with cluster algebras. The appropriate reduction of the Lax pair 
(j46p for discrete KdV, via the staircase method, was the key to finding explicit expressions for first 
integrals, and two ways were presented to prove that these are in involution. The corresponding reduc- 
tion of the lattice equation ([3]) was also seen to be completely integrable (and even superintegrable), 
with additional first integrals appearing due to the presence of the d-integral J„. 
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An interesting feature of all these reductions is that, although they are autonomous difference 
equations, they have various difference equations with periodic coefficients associated with them, such 
as (IMD, (IMl) and (l671) . 

There are several ways in which the results in this paper could be developed further. In particular, 
it would be interesting to understand the Poisson brackets for the reductions in terms of an appropriate 
r-matrix. It would also be instructive to make use of the bi-Hamiltonian structure to perform separation 
of variables, by the method in [3], and this could be further used to obtain the exact solution of these 
difference equations in terms of theta functions associated with the hyperelliptic spectral curve ([5T]) . 
Finally, it would be worthwhile to extend the results here to (si, — ■S2)-reductions of discrete KdV, and 
to reductions of other integrable lattice equations that have not been considered before. 
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Appendix: Proof of Theorem [T] 



In order to prove formula (jl4p . we note that 

{Un-d,Un-r-l} = {[qi;Pl " + " —,P2 + 92 + 93, • • • ,Pd--l + Qd-l + qd,Pd], 

Pd 

[qd-r;Pd-r + qd~r + qd^r+l,Pd~r+l + Q'd-r+l + qd-r+2, ■ ■ ■ ,Pd-l + qd-1 + QdiPd]} 

E, dUn-d dUn^r-1 _ dUn-d dUn-r-l -. 
dpi dqi 



i>d—r 

Then, for i > d — r > 2, we have 



dqi dpi 



du 



n—d 



dpi 
dpd 

dUn-d 

dqi 

dUn-d 

dqd 

dUn-r-l 

dpi 

dUn-r-l 

dqi 



-iy{Un-d - Unf ■ ■ ■ {Un-d+i-1 " Un+i-lf , i < d, 

-iy{Un-d - Unf ■ ■ ■ {Un-1 - Un+d-lf " {Un~d - Unf {Un-l - Un+d-lf 
-iy^^{Un-d - Unf ■ ■ ■ {Un-d+i-2 - M„+i_2)^(l - {Un-d+i-l - Un+i-lf) , i < d 
Un-d - Unf + {-lY~^{Un-d " Unf ■ ■ ■ {Un-2 " Un+d-2f , 
_lY+^ 1 '^{un-T-l — Un+d-r-iy ■ ■ ■ (Un-d+i-l — Un+i-l)^ , 

— 'i-y'^^~'^{Un-r-l — Un+d-r-iy ■ ■ ■ {Un-d+i-2 — 'Un+j-2)^(l — {Un-d+i-1 — ^tn+j-l)^), 



where the last identity holds for i < d, and 



dUn~r~l f -\\T/ \2 / \2 
= l~lj [Un-r-l — Un+d-r-l) ■ ■ ■ {Un-2 — Un+d-2) , 

dqd 



du 



n—d 



from which we obtain {un-d^Un-r-i} = 0. Next, we use 

{Un-d,Un+j} = {Un-d,qj+l} o 

CPj+1 

where < j < ci — 1, from which the formulae (jlSp and (jl6p follow. 
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